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ABSTRACT 

SAR  is  widely  used  today  in  earth  observation  and  is  a  potential  means  for  military  surveillance.  However, 
systems  only  operating  in  the  basic  SAR  mode  will  be  superseded  by  new  developments  which  fulfil  much 
more  user  demands.  Highest  resolution  in  flight  direction  combined  with  a  far  range  operation  cannot  be 
reached  in  stripmap  mode.  The  spotlight  mode  has  to  be  implemented  demanding  for  special  capabilities 
of  the  hardware  and  the  data  processing  chain.  To  achieve  a  fine  resolution  in  range  direction,  transmit 
signals  with  very  large  bandwidth  have  to  be  used.  Technological  imitations  force  the  developers  to 
implement  new  strategies  to  acquire  the  necessary  bandwidth. 

Multichannel  SAR  systems  offer  great  advantages  compared  to  the  conventional  ones.  Depending  on  the 
antenna  configuration  such  systems  allow  to  measure  the  height  of  image  pixels  or  to  detect  moving 
objects  in  the  illuminated  scene. 


1.0  VERY  HIGH  RESOLUTION  SAR 


1.1  Relation  of  resolution  and  k-set 


The  analyses  of  the  resolution  of  a  SAR  system  must  distinguish  range  and  azimuth  (flight  direction) 
resolution.  While  a  broad  bandwidth  of  the  transmitted  and  received  pulses  is  responsible  for  a  fine  range 
resolution,  the  azimuth  resolution  can  be  improved  by  using  a  short  wavelength  or  by  a  large  processing 
angle.  However,  in  both  directions  the  designated  measures  cause  similar  effects  in  the  wave  number 
domain  representation  of  the  echo  data.  Thus,  a  common  examination  of  the  relation  between  resolution 
and  k-set  can  be  accomplished. 
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Figure  1:  Relation  of  k-set  and  image  resolution.  Left  diagram:  representation  of  signals  in  spatial  and  wave 
number  domain.  Right  diagram:  InUuence  of  the  point  spread  function. 
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Figure  1  illustrates  the  relation  of  signals  in  the  spatial  domain  and  its  Fourier  correspondences  in  k-  or 
wave  number  domain.  The  left  diagram  depicts  the  complete  one-dimensional  image  formation  process: 

•  a  given  scene  will  be  illuminated  using  an  arbitrary  transmit  pulse,  usually  a  chirp  signal 

•  reflected  echoes  of  all  scene  targets,  two  chirp  signals  in  the  example,  will  be  received 

•  a  matched  filter  is  used  for  compression  by  a  factor  determined  by  the  time-bandwidth  product  of 
the  transmit  signal  (signal-to-noise-ratio  is  improved  accordingly). 

The  reconstructed  image  is  the  convolution  output  of  the  scene  and  the  point  spread  function.  The  right 
diagrams  represent  a  generalized  view  of  the  same  relationship.  The  image  properties  are  completely 
defined  by  the  reflectivity  and  the  point  spread  function.  The  concrete  appearance  of  the  transmit  signal  is 
irrelevant  for  the  imaging  process.  A  corresponding  situation  results  in  azimuth  direction.  A  k-set  will  be 
generated  by  observation  of  the  scene  from  different  directions.  A  large  observation  angle  thus  produces  a 
large  k-set  and  a  fine  azimuth  resolution.  Since  the  point  spread  function  is  the  Fourier  correspondence  of 
the  indicator  function  in  k-space,  it  is  clear,  that  the  wave  number  domain  has  to  be  covered  without  gaps. 
Otherwise,  high  sidelobes  arise  in  the  target  response.  Figure  2  illustrates  the  point  spread  function  for 
different  cases  of  k-sets.  Case  a)  is  typical  for  a  SAR  with  a  limited  bandwidth.  The  observation  angle  is 
also  bounded,  perhaps  by  operation  in  stripmap  mode.  K-set  acquisition  in  case  b)  is  done  using  a 
monofrequent  radar  system  and  an  overall  observation.  Though  the  bandwidth  of  the  transmit  signal  is 
negligible,  resolution  can  be  obtained  in  all  directions.  Sidelobes  increase  because  of  the  small  proportion 
of  the  covered  k-area.  Case  c)  depicts  the  situation  of  a  system  with  very  large  relative  bandwidth  and  a 
wide  processing  angle.  This  configuration  produces  problems  in  the  data  processing,  because  very 
accurate  algorithms  have  to  be  used.  Diagram  d)  shows  the  effect  when  undersampling  occurs  in  the  wave 
number  domain  which  results  in  sidelobes.  Moreover,  the  distance  of  samples  in  k-space  is  reciprocal  to 
the  image  size.  If  the  illuminated  scene  (according  to  the  antenna  beam  angle)  is  larger  than  this  image 
area,  aliasing  appears.  Compared  to  this  case,  the  unconformable  structure  of  the  k-set  in  diagram  e)  leads 
to  reduced  sidelobes.  However,  the  integrated  sidelobe  level  remains  high.  Aliasing  can  be  avoided  by  this 
measure.  Diagram  f)  shows  a  well  formed  k-set  with  a  superior  point  spread  function. 

1.2  Measures  to  enhance  range  resolution 
1.2.1  Overview  of  methods 

With  the  increasing  demands  according  to  the  resolution  capability  of  modem  SAR  systems,  the  range 
resolution  encounters  limitations.  The  necessary  bandwidth  claims  for  superior  Antenna,  HF-  and  central 
electronics  and  A/D  converter. 

Often,  the  A/D  conversion  of  the  received  echoes  plays  a  limiting  role  which  also  influences  the  data 
processing.  For  a  resolution  of  a  decimetre,  a  bandwidth  of  1.5  GHz  has  to  be  used.  However,  even 
modem  A/D  converters  with  adequate  quantization  offer  sampling  rates  of  only  800  MHz.  Consequently, 
a  direct  sampling  of  the  echoes  cannot  be  realized.  Two  methods  allow  to  acquire  data  with  the  demanded 
bandwidth. 

The  deramping  is  implemented  in  a  couple  of  current  radar  systems  to  reduce  the  signal  bandwidth  which 
has  to  be  recorded.  The  application  is  not  restricted  to  SAR,  but  also  ISAR  systems  like  TIRA  (tracking 
and  imaging  radar)  of  FGAN  apply  this  method.  The  idea  is  to  transmit  linear  frequency  modulated  pulses 
and  to  convert  the  received  echoes  using  an  equivalent  modulated  reference  signal  to  the  base  band.  If  the 
frequency  modulation  of  the  reference  signal  compensates  for  the  echo  modulation,  the  bandwidth  of  the 
acquired  signal  will  be  reduced  considerably. 
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Figure  2:Examples  of  different  shaped  k-sets  and  their  corresponding  point  spread  functions,  a)  small  angular 
sector,  small  bandwidth,  b)  Large  angular  sector,  small  bandwidth,  c)  Large  angular  sector,  large  bandwidth, 
d)  undersampled  observation  angle,  e)  scattered  frequency  band,  f)  well  formed  k-set,  tuned  frequency  bands 

The  alternative  solution  is  to  share  the  total  bandwidth  between  some  pulses  having  different  centre 
frequencies.  These  pulses  will  be  transmitted  and  received  subsequently  and  can  be  digitized  by  the 
available  sample  rate.  In  the  processing  process,  the  partial  signals  will  be  combined  coherently. 

1.2.2  Deramping 

The  idea  of  deramping  is  to  reduce  the  signal  bandwidth  within  the  analogue  section  of  the  receiver  and  to 
perform  the  digital  signal  processing  with  a  low  sampling  rate.  If  linear  frequency  modulated  transmit 
signals  will  be  used  the  down  conversion  with  an  also  linear  FM  reference  signal  leads  to  mono-frequent 
signal  contributions.  Figure  3  illustrates  this  subject. 
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Figure  3:  Illustration  of  the  dramping  method:  The  upper  diagram  represents  the  reflectivity  of  a  scene 
having  three  point  scatterers.  Beneath,  the  transmit  signal  and  scene  echoes  are  depicted.  The  lower  diagram 
shows  the  received  signal  after  deramping.  Now,  echoes  are  mono  frequent  contributions  of  the  total  signal. 

The  indicated  low  pass  filter  causes  a  range  limitation. 

Some  characteristics  become  clear  by  this  presentation: 

•  The  reference  signal  must  be  a  scene  size  larger  than  the  transmit  pules. 

•  The  larger  the  imaged  range  area  should  be,  the  larger  the  bandwidth  of  the  reference  signal  has  to 
be  chosen. 

•  After  deramping,  the  total  bandwidth  is  proportional  to  the  scene  extension.  A  low  pass  filtering  of 
the  received  signal  corresponds  to  a  limitation  of  the  received  range  interval.  The  signal 
bandwidth  thus  can  be  adopted  to  the  available  samping  rate. 

Now,  we  will  analyse  the  deramping  mathematically.  The  linear  frequency  modulated  transmit  signal 

st(r)  =  rect  ei(fcr'or_^r2) 

is  assumed  to  be  given  which  covers  the  interval  [kr  0~Kr  /2,kr0+Kr/2\.  Kr  is  the  size  of  the  wave 

number  interval  and  y  denotes  the  wave  number  slope  in  r  .  The  length  of  the  signal  is  Kr  /  y .  If  this 
signal  will  be  reflected  by  a  point  target  with  reflectivity  a  at  distance  r  ,  the  sensor  receives  the  echo 

z(r)  =  st(r  —  r')  =  a  rect  ^  ^  J  ^  ej[kr,o(r~r  )-i(r_r  )2] 

The  deramping  will  be  performed  by  down-conversion  using  the  reference  signal 

znf(r)  =  e»[^°<r-r‘>)-Fr-r«)a] 
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where  rc  denotes  a  reference  range  within  the  range  of  the  image.  The  deramped  signal  results  in 


z(r)zre/  (r) 


=  a  rect 


r  —  r 


Kr/  7 


e-jkr,o(r,-rc)e-j%(r'2-i2)ej'y(r,-rc)r 


The  last  exponential  term  possesses  a  linear  phase  in  r  which  depends  on  the  location  r'of  the  imaged 
point  target.  After  a  Fourier  transform  this  term  causes  a  shift  in  the  wave  number  domain.  It  follows 


Zd  (J^r) 


aKr 


e-j(fcr,0+fcr)(r'-rc)e-jfcrrCej§(r'-rc)2si  J 


»  {^(Av_7(r,_rc))}- 


The  first  exponential  term  supports  the  accurate  information  of  the  range  difference  rf  -  r  ,  which  is  the 

basis  for  the  coherent  processing  of  pulses  in  flight  direction.  The  second  phase  term  is  independent  of  the 
object  distance  and  can  be  compensated  easily.  The  third  term,  known  in  the  literature  as  residual  video 
phase,  has  to  be  examined  and  removed,  if  necessary.  An  exact  solution  is  given  in  [2]  and  an 
approximation  is  presented  in  [15].  Finally,  the  si-function  represents  the  amplitude  and  shows  that  the 
signal  has  been  compressed.  The  resolution  can  be  determined  to  Sr  =  2n  /  Kr  which  is  the  common 
relation  of  resolution  and  the  size  of  the  k-set. 


1.2.3  Synthetic  bandwidth 

Another  idea  to  overcome  the  limitation  of  the  bandwidth  is  to  transmit  a  series  of  signals  which  cover 
different  frequency  intervals.  The  coherent  combination  of  the  received  echoes  leads  to  the  final 
resolution.  The  principle  of  synthetic  bandwidth  (also  denoted  as  stepped  chirp)  is  depicted  in  Figure  4. 
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Figure  4:  Principle  of  synthetic  bandwidth:  The  upper  images  show  three  transmit  signals  covering  different 
frequency  intervals  and  their  echoes.  The  second  row  depicts  the  spectra.  The  last  diagram  shows  the 
combined  signal  covering  the  whole  wave  number  interval. 
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The  mathematical  analyses  starts  assuming  transmit  signals  stn(r )  for  n  =  1, . . .n  .  If  st  n  (r)  has  constant 

power  density  L  in  \kr  0  -  nKr  /  2  + (n  —  l)Ar ,  kr  0  -  nKr  /  2  +  nKr  J  received  echoes  of  a  scene  with 
reflectivity  a{r )  after  matched  filtering  result  as 


^771,71  (^r )  L  lcct 

Coherent  superposition  leads  to 

ft 

Zm,Y,  (kr)  =  ^  ^  Zm,n  (kr)  =  A  Tect 

71—1 


—  [^r,o  —  (n  —  l)A"r/2  +  (n  —  l)ATr] 


a  (Av). 


kr  kro\  • 


The  inverse  Fourier  transformation  into  the  spatial  domain  results  in 

*m, E(0  =  (PE  *  a)(0> 


where 


Ps(r)  = 


n/C 


is  the  point  spread  function.  Compared  to  the  resolution  of  one  single  pulse  (with  partial  bandwidth)  the 
resolution  will  be  S  =  2tc  l{nKr ) ,  which  is  enhanced  by  a  factor  of  n  . 


2.0  INTERFEROMETRY 


2.1  Two-antenna  interferometry 

In  a  conventional  SAR  image  a  scene  point  will  be  depicted  according  to  its  distance  to  the  flight  path  and 
the  location  in  flight  direction  where  the  point  is  perpendicular  to  the  synthetic  aperture.  Thus,  reflecting 
points  at  circles  around  the  flight  path  will  be  located  in  the  result  in  one  single  image  pixel.  It  is  not 
possible  to  localise  the  point  in  three  dimensions.  In  terrain  height  retrieving  interferometric  SAR  (IfSAR) 
the  elevation  angle  of  each  ground  pixel  is  estimated  by  determining  the  phase  difference  of  the  signals 
from  at  least  two  antennas  separated  along  the  cross-track  direction. 


Zy=ae"^rf 


a|VjP(r>_r2) 

(3:  Wavenumber 


Figure  5:  Geometry  of  single-pass  cross-track  SAR  interferometry. 

Figure  5  depicts  the  geometry  for  a  two  antenna  configuration.  Knowing  the  ground  pixel  and  the 
antennas'  positions,  the  terrain  height  is  obtained  from  the  elevation  angle.  Due  to  the  2k  ambiguity  of 
this  phase  difference,  a  reference  point  is  needed  to  determine  the  phase's  principal  value  and  the  allowed 
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terrain  height  variation  is  strictly  limited.  This  limitation  can  be  avoided  by  a  suited  phase  unwrapping 
procedure  (i.e.  adding  the  appropriate  number  of  2 n  to  each  pixel  to  minimize  discontinuities),  assuming 
the  terrain  slope  not  to  exceed  a  certain  extent,  in  order  to  avoid  phase  aliasing.  Due  to  decorrelation 
effects  this  phase  unwrapping  is  an  error-prone  process,  leading  partly  to  poor  height  estimations.  So, 
phase  unwrapping  is  the  key  component  of  SAR  interferometry  processing,  if  no  digital  elevation  model  is 
used.  Most  phase  unwrapping  algorithms  can  be  categorised  into  path-following  and  least-squares 
algorithms. 

An  analyses  of  the  k-set  of  acquired  interferometric  data  shows  that  a  third  dimension  has  also  been 
established  in  the  wave  number  domain.  Figure  6  illustrates  this  fact.  The  Fourier  counterpart  thus 
incorporates  resolution  capabilities  in  all  directions. 


Figure  6:  Constellation  of  interferometric  SAR  acquisition.  The  pair  of  k-sets  of  the  two  antennae  captures  the 

third  dimension  in  k-space. 

A  further  source  for  phase  corruption  is  given  by  contributions  of  different  elevation  angles  within  one 
resolution  cell  (so-called  layover).  Typical  problems  of  that  kind  rise  up  in  many  images  with  bridges  over 
a  valley  or  a  similar  constellation.  In  such  a  case,  the  reflectivity  of  multiple  scatterers  superpose  each 
other  and  in  the  image  pixel  only  the  sum  of  all  appears.  Hence,  the  phase  which  carries  the 
interferometric  height  information  of  the  pixel  is  a  mixture  of  the  phase  of  all  reflecting  elements  sharing 
the  same  image  pixel.  If  one  contribution  dominates  the  others,  its  height  information  will  prevail. 
Otherwise,  the  calculated  height  will  be  useless. 

Figure  7  depicts  an  interferometric  SAR  image  which  has  been  recorded  using  two  antennae. 

2.2  Multi  baseline  interferometry 
2.2.1  Principle 

The  height  estimation  accuracy  is  inversely  proportional  to  the  baseline  of  the  interferometer.  But, 
increasing  the  baseline  results  in  higher  fringe  frequencies,  thus  making  the  phase  unwrapping  still  harder. 
Therefore,  it  is  quite  natural  to  reduce  these  problems  by  using  a  multi  baseline  system:  the  ambiguities 
can  be  removed  by  a  suitable  small  baseline,  thus  retrieving  the  absolute  phase  with  no  more  need  for 
reference  points  and  phase  unwrapping,  and  the  height  estimation  accuracy  can  be  sustained  by  a  longer 
baseline.  As  a  second  advantage,  multi  baseline  IfSAR  allows  for  estimation  of  the  different  elevation 
angles,  caused  by  mixing  different  height  pixels  within  one  resolution  cell,  by  means  of  superresolution 
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methods. 

In  the  literature,  some  theoretical  investigations  or  simulations  of  the  multi  phase  centre  interferometry  can 
be  found  [3,9,10]. 


Figure  7:  Perspective  view  of  an  interferometric  SAR  image. 

Following  on,  we  analyse  an  antenna  configuration  with  three  single  antennae  separated  in  cross-track 
direction.  One  could  expect  that  a  very  short  baseline  results  in  ambiguity  and  phase  unwrapping  removal. 
On  the  other  hand,  we  have  to  minimize  the  sidelobe  level  of  the  three-antenna  array,  because  high 
sidelobe  levels  could  result  in  wrong  elevation  estimates  when  using  noisy  data.  This  minimisation 
determines  the  position  of  the  third  antenna.  Figure  8  depicts  the  scan  pattern  of  a  three-antenna 
configuration  with  baseline  ratios  1:2:3  compared  to  a  two-antenna  configuration.  As  can  be  seen  the 
ambituity  interval  is  tripled  compared  to  the  two-antenna  array.  An  irrational  ratio  would  in  fact  result  in 
an  infinite  ambiguity  interval  at  the  cost  of  increasing  sidelobe  levels.  However,  this  could  only  be 
exploited  with  noise-free  data. 


Scan  Pattern 


A  1  "  1  "  0  "  1  ^  '  5 
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Figure  8:  Scan  pattern  of  two-antenna  (dotted)  and  three-antenna  (solid)  configuration  with  baseline  ratios 

1:2:3 

One  possibility  to  exploit  the  additional  information  of  the  third  antenna  is  to  support  the  conventional 
phase  unwrapping  processes  of  one  antenna  pair  by  the  additional  information  of  the  third  antenna. 
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However,  a  vector  based  interpretation  of  the  three  measured  complex  valued  image  pixels  of  each 
resolution  cell  provides  better  results  and  additional  features. 

2.2.2  Pixel  based  elevation  estimation 

We  estimate  the  elevation  angle  of  a  pixel  in  principle  independent  of  any  neighbour  pixels  on  the  basis  of 
a  combined  evaluation  of  the  three  measured  signals.  The  following  model  can  be  used  for  a  pixel  v  of 
the  image,  which  has  three  components  according  to  the  tree  antennae. 

Zv  =Cvav(Uv)+nv 

Here,  c  is  the  complex  SAR  image,  u  is  a  unit  vector  to  the  analysed  pixel,  n  is  uncorrelated  receiver 
noise  and 


1 

Vs, 

_ 1 

a(u)  = 

/^2(u)e//?l,"'^, 

L(u)ejMa^ 

where  r.  denotes  the  position  of  sensor  i  and  kt  comprises  the  two-way  characteristics  as  well  as  the 
relative  phase  shifts  of  the  i  -th  transmit/receive  pair.  a(u)  is  called  DOA-vector  (direction  of  arrival). 

In  a  first  step,  a  calibration  of  the  system  has  to  be  carried  our,  i.e.  k(u)  has  to  be  determined.  After  flat 
earth  compensation  (i.e.  removing  a  slant  range  dependent  phase  to  get  the  expected  phase  term  for  a 
topographically  flat  earth  surface),  an  estimation  for  k  is  determined  by  the  eigenvector  corresponding  to 
the  largest  eigenvalue  of  the  covariance  matrix  given  by  the  flat  earth  compensated  signal  vector  zv .  This 

leads  to  a  calibrated  representation  of  the  signal,  which,  after  some  averaging,  serves  as  input  to  a 
monopulse  estimator.  The  thus  estimated  Au  determines  the  deviation  of  the  flat  earth  and  allows  to 
calculate  the  height  pixel-by-pixel  by  simple  geometrical  operations. 

If  we  look  at  the  eigenvalue  distribution  of  the  two  largest  eigenvalues  \  and  X2  of  the  averaged  signal, 
we  expect,  with  the  exception  of  some  leakage  out  of  the  signal  subspace  by  high  intensities,  a  high 
\  /  A2  ratio,  if  only  one  single  Au  contributes  to  the  pixel.  On  the  other  hand,  a  moderate  ratio  will  result 

if  two  Au  ,  belonging  to  different  heights,  are  mixed  up  within  the  same  pixel.  Thus,  we  are  able  to  decide 
weather  one  pixel  contains  contributions  from  more  than  one  scene  height.  The  eigenvectors  can  be  used 
in  a  subsequent  step  to  estimate  the  heights  which  belong  to  the  current  image  pixel. 


3.0  MOVING  TARGET  INDICATION  (MTI) 

3.1  Properties  of  moving  target  echoes 

Moving  objects  pose  problems  in  conventional  SAR  since  the  relative  motion  between  the  radar  sensor 
and  the  scene  behaves  no  longer  as  presumed  originally.  SAR  processing  of  acquired  moving  target 
echoes  leads  to  defective  images  of  the  targets:  misplaced  in  the  resulting  image,  defocused  or  not 
appearing  at  all. 

The  problem  can  be  formalised  similar  to  SAR  imaging.  A  sensor  moves  along  the  x  -axes  with  a  constant 
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velocity.  A  target  at  position  {xt,yt )  moves  in  arbitrary  direction  with  also  constant  velocity.  The  distance 

and  the  direction  from  the  sensor  to  the  target  depend  on  slow-time  and  the  target  parameters  i9  .  Figure  9 
depict  the  situation. 


Figure  9:  Geometry  for  a  SAR  constellation  with  a  moving  target. 

Analysing  recorded  raw  data  in  slow-time  incorporates  different  parameters:  time,  range  and  frequency. 
While  two-dimensional  presentations  of  the  data  in  range  and  frequency,  range  and  time  as  well  as  in  time 
and  frequency  are  commonly  used,  a  three  dimensional  illustration  like  depicted  in  Figure  10  is  seldom 
presented.  The  appearance  of  the  relation  depends  on  all  target  parameters.  However,  the  parameters 
cannot  be  unambiguously  determined  from  the  course  of  the  data. 


Figure  10:  Three-dimensional  SAR  data  representation. 

Consequences  of  unconsidered  motion  of  targets  in  conventional  SAR  imaging  is  illustrated  in  Figure  11. 
Radial  motion  only  causes  a  shift  of  the  Doppler  band.  The  slope  of  the  Doppler  history  remains 
unchanged,  thus  the  matched  filter  compresses  the  echo  of  the  target  to  a  small  object.  Because  of  the 
shifted  Doppler  spectrum  a  portion  of  the  signal  energy  will  not  contribute  to  the  image  so  the  target 
response  will  be  decreased.  Since  the  Doppler  modulation  is  approximately  linear,  the  shift  in  Doppler 
domain  corresponds  to  a  shift  in  range  direction.  The  radial  motion  thus  also  results  in  a  range  offset  of  the 
imaged  target.  Echoes  of  a  target  which  moves  parallel  to  the  flight  path  have  a  Doppler  slope  which 
deviates  from  the  slope  of  clutter  objects.  The  commonly  used  matched  filter  is  not  able  to  compress  such 
a  target,  because  uncompensated  phase  fluctuation  remains  after  the  matched  filter  process.  Anti-parallel 
motion  leads  to  a  similar  effect.  If  a  target  has  a  radial  component  which  shifts  the  Doppler  frequency 
history  out  of  the  processed  clutter  band,  it  will  not  be  imaged. 
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Figure  11:  Effects  of  motion  in  SAR  imaging. 

A  couple  of  heuristic  approaches  exist  to  cope  with  moving  targets  using  a  conventional  single  channel 
SAR  system: 

•  Detection  outside  the  clutter  bandwidth 

•  Chirp  slope  analysis 

•  Chirp  filter  bank 

•  MTI  via  multilook  images  (change  detection) 

•  Range  walk  detection 

•  Chirp  offset  detection 

•  Reflectivity  displacement  method 

•  Relation  to  context: 

•  Shadow  detection 

•  Off-road  detection 

However,  these  methods  exhibit  weaknesses  particularly  for  fast  and  slowly  moving  targets  and  for  low 
SNR.  Efficient  MTI  capability  is  only  obtainable  using  multiple  channels  in  flight  direction. 

3.2  Multichannel  MTI 

Figure  12  illustrates  the  idea  of  multichannel  MTI  and  the  realisation  concept  in  the  SAR  context.  Image 
a)  depicts  the  displaced  phase  centre  antenna  technique  (DPCA):  two  antennae  installed  at  different 
positions  of  the  sensor  platform  acquire  radar  data,  using  the  same  geometry  (determined  by  the  antenna 
position  at  pulse  acquisition  instant)  but  at  different  moments.  Obviously,  fixed  targets  produce  identical 
echoes,  thus  they  cancel  out  in  the  difference  signal.  Remaining  energy  is  caused  by  objects  with  radial 
motion  and  noise.  DPCA  is  very  illustrative  but  can  only  be  used  in  extended  variants.  SAR  data  of  a  two- 
antenna  system,  like  it  is  depicted  in  b),  can  be  exploited  by  along-track  interferometry.  In  a  first  step,  a 
complex  SAR  image  will  be  calculated  for  each  channel.  Afterwards,  the  images  will  be  multiplied 
complex  conjugated.  The  remaining  phase  indicates  radial  motion  of  moving  targets.  For  fixed  targets  the 
resulting  phase  vanishes.  A  result  of  along-track  interferometry  is  shown  in  Figure  13.  To  improve  the 
MTI  capability,  more  antennae  have  to  be  used.  Figure  c)  depicts  this  technique.  The  acquired  data 
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represent  samples  of  the  wave  field  in  space  and  time.  With  the  aid  of  space  time  adaptive  processing 
(STAP)  a  couple  of  motion  information  can  be  calculated  out  of  the  data. 
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Figure  12:  Ideas  for  multichannel  MTI  processing:  a)  displaced  phase  centre  antenna  (DPCA),  b)  application 
of  DPCA  to  SAR  case,  c)  extension  of  the  principle  to  space  time  adaptive  processing  (STAP) 


Figure  13:  Along-track  interferometry  of  the  ship  lock  Brunsbiittel.  The  image  colour  is  chosen  according  to 

the  differential  phase  of  the  two  SAR  images. 
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